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There has been a rapidly growing interest on the interplay between spin-orbit coupling (SOC) and Hubbard 
interaction U in correlated materials. A current consensus is that the stronger the SOC, the smaller is the critical 
interaction U c required for a spin-orbit Mott insulator, because the atomic SOC splits a band into different total 
angular momentum bands narrowing the effective bandwidth. It was further claimed that at large enough SOC, 
the stronger the SOC, the weaker the U c because in general the effective SOC is enhanced with increasing 
electron-electron interaction strength. Contrary to this expectation, we find that, in orthorhombic perovskite 
oxides (Pbnm), the stronger the SOC, the bigger the U c . This is originated from a line of Dirac node in J e ff =1/2 
bands near the Fermi level inherited from a combination of the lattice structure and a large SOC. Due to this 
protected line of nodes, there are small hole and electron pockets in SrlrOj, and such a small density of states 
makes Hubbard interaction less efficient in building a magnetic insulator. The full phase diagram in U vs. 
SOC is obtained, where non-magnetic semimetal, magnetic metal, and magnetic insulator are found. Magnetic 
ordering patterns beyond U c are also presented. We further discuss implications of our finding in relation to 
other perovskites such as SrRhOi and SrRu0 3 . 

PACS numbers: 



INTRODUCTION 

Perovskite oxides with the chemical formula AMO3 where 
A is a cation and M is a transition metal, exhibit an exception- 
ally wide range of properties including anomalous Hall effect, 
colossal magnetoresistance, ferroelectricity, ferromagnetism, 
and superconductivity. Such an ample variety in a rather sim- 
ple structure indicates that a detailed balance between charge, 
spin, structure, and correlation is important in determining the 
ground state. 

In particular, orthorhombic perovskite (point group symme- 
try, Pbnm) oxides are a large class of anisotropic oxides based 
on AMO3 where MO;, octahedra are distorted from the sym- 
metric cubic structure. Among them, SrRuC>3, SrRhC>3 and 
SrIrC>3 (called perovskite ruthenates, rhodates, and iridates re- 
spectively), display correlated metallic ground states. How- 
ever, their magnetic properties differ hinting a crucial role of 
electron interaction. SrRuC>3 is a ferromagnetic metal |Q}[3) 
and SrRhC>3 a metal near a critical point IUJ6), while Srlr03 
is a semimetal with an extremely small number of charge car- 
riers without any magnetic moment [7-9 1. Given that Ir has 5d 
orbitals in the outer shell, while Rh and Ru have 4d orbitals, 
Hubbard interaction is expected to be smaller in iridates IflOl . 
Indeed it was found that the optical gap due to Hubbard in- 
teraction is about 0.5eV in Sr?IrO,i lfrn . a sister compound of 
SrIrC>3. This leads to a naive conclusion that iridates should 
be better metal than rhodates or ruthenates, but the reality is 
the opposite. 

What is missing in the above discussion is the SOC. Ir is 
heavier than Rh and the SOC strength is comparable to Hub- 
bard interaction in iridates ifTTTl . Since the atomic SOC is a 
local interaction, the electronic energy level splits into differ- 



ent total angular momentum J levels. For example, starting 
from the atomic limit, five d-orbitals split into and e CJ lev- 
els due to the octahedral crystal field, and t2g further splits 
into J e ff = 3/2 and J e ff = 1/2 via the SOC when the crys- 
tal field splitting is larger than the strength of SOC. Once 
these bands form, a larger SOC leads to a smaller bandwidth 
of J e ff =1/2 separated from J e ff = 3/2. Thus, the larger 
the SOC, the larger the ratio between Hubbard interaction^) 
and the bandwidth(W), U/W where W is the bandwidth of 
Jeff = 1/2. While the absolute strength of U is smaller in 
iridates, its effect (given by the ratio U/W) is amplified. This 
is indeed observed in a layered perovskite, Sr2lr04, dubbed 
a spin-orbit Mott insulator lfl2llT8ll . To explain the metallic - 
ity of Srlr03 compared to insulating Sr2lr04, it was further 
suggested that Srlr03 has a larger bandwidth comparing to 
quasi-two dimensional Si"2lr04 |[T9l |20| . A growing consen- 
sus is that the larger the SOC, the smaller the critical interac- 
tion strength U c that is required for the phase transition from 
metal to Mott insulator lfP3ll2TI 

However, once the SOC splits the bands into different 
J e ff bands, its effect on the bandwidth of J e ff — 1/2 is min- 
imal, and the interplay between the SOC and the electron- 
electron interaction is intriguing. It was claimed that in gen- 
eral the effective spin-orbit coupling is enhanced with in- 
creasing strength of the electron-electron interaction leading 
to the same conclusion that a larger SOC leads to a smaller 

In this paper, we show a counter example where the com- 
mon wisdom does not apply. We study the interplay between 
SOC and Hubbard interaction in orthorhombic perovskite ox- 
ides (Pbnm). It is found that the bigger the SOC, the larger 
the U c in orthorhombic perovskites inherited to the lattice 
structure. When SOC is moderate (close to the true SOC 
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FIG. 1: Crystal structure of orthorhombic perovskite SrlrOv Sr, Ir 
and O atoms are shown in aqua, blue and red. The octahedra shown 
are rotated about the z-axis and tilted about [110]-axis making the 
unit cell four time bigger than that of the cubic perovskite structure. 

in SrIrOs), the band dispersion exhibit a line of Dirac node 
protected by the symmetry of the lattice. We propose that 
semimetallicity in SiirC>3 compared to insulating Si"2lr04 is 
due to such a small density of states, which in turn requires a 
larger U c for the transition to a Mott insulator. Hubbard in- 
teraction in iridates is smaller than this U c , and thus SrIrC>3 
remains metallic with small Fermi pockets . Beyond U c , non- 
collinear and non-coplanar magnetic structures appear, and 
the overall phase diagram contains ferromagnetic metal, non- 
magnetic semimetal and magnetic insulator. Below we will 
show the band structures computed for SrIrC>3, where we use 
Hubbard U and SOC strength a as tuning parameters to un- 
derstand different phases realised in other orthorhombic per- 
ovskite oxides such as SrRuC>3 and SrRh03. Our findings 
suggest that the SOC together with Hubbard interaction U 
plays an important role in realising different ground states in 
SrRu0 3 I22H27I. SrRh0 3 @, and Srlr0 3 J28H30|. 

The paper is organised as follows. In the following section, 
the details about the crystal structure is presented. In Sec. 3, 
computational method is explained, and the band structures 
and phase diagram in U vs. SOC are presented in Sec. 4. 
Magnetic metal and insulator appear at small and large U re- 
spectively, and their magnetic ordering patterns depend on the 
SOC which will be shown in Sec. 5. A brief summary and 
implications of our findings are listed in the final section. 



CRYSTAL STRUCTURE 

Fig. [T] shows the crystal structure of the orthorhombic per- 
ovskite Srlr03 with Sr, Ir and O atoms as aqua, blue and red 
balls. As can be seen in Fig.[T] the octahedra enclosing the 
Ir atoms are rotated about the z-axis and tilted about [110]- 
axis. Due to these rotations and tilts, there are four formula 



units of Srlr03 in a unit cell and the octahedra also get dis- 
torted. For any two connected octahedra, the rotations are in 
the same (opposite) direction if the two enclosed Ir atoms lie 
in different (the same) layers, whereas the tilts are opposite for 
all nearest neighbour octahedra. 

The experimental lattice parameters of this Pbnm phase 
of Srlr0 3 are a = 10.5136 a.u., b = 10.5688 a.u. 
and c = 14.9 a.u., and an asymmetrical unit is: a Sr 
at (0.5085,0.4901,0.25), an Ir at (0.5,0,0) and two O at 
(0.506, 0.073, 0.25) and (0.292, 0.714, 0.044) OTI . 

This structure is primitive orthorhombic for which the sym- 
metry elements include two b glide planes perpendicular to 
x-axis at x/a = 1/4 & 3/4, two n glide planes perpendicular 
to j/-axis at y/b — 1/4 & 3/4 and two mirror planes perpen- 
dicular to z-axis at z/c = 1/4 & 3/4. Here, a b (n) glide 
plane means that a reflection across the plane followed by a 
translation of a/2 ([a + c]/2, i.e, along the diagonal) trans- 
forms the structure to self coincidence. Furthermore, there 
are four 2\ screw axes parallel to each of the three primi- 
tive lattice vectors a,b and c. The 2i screw axes parallel to 
a or x-axis are at (y/b, z/c) = (1/4, 0), (1/4, 1/2), (3/4,0) & 
(3/4,1/2); those parallel to b or j/-axis are at (x/a, z/c) = 
(1/4, 1/4), (1/4, 3/4), (3/4, 1/4) & (3/4,3/4); and those par- 
allel to c or z-axis are at (x/a,y/b) = (0,0),(0, 1/2), (1/2,0) 
& (1/2, 1/2). 

There are eight inversion centres at x/a, y/b, z/c e {0, 1/2). 
The four Ir atoms in the unit cell sit at four of these. This also 
means that all the octahedra in Fig.[T]are inversion symmetric. 
While this is obvious in case of a cubic perovskite structure 
which forms regular octahedra around the Ir atoms, it is not 
so in this case where the octahedra are distorted. Two of the 
2i screw axes parallel to c passes through the Ir atoms. These 
screw axes and the four inversion centres at Ir locations are 
necessary for the existence of the mirror planes at z/c = 1/4 
& 3/4, which connect the octahedra in two different layers 
through the reflection symmetry. It was found in Ref. 11321 
that breaking this mirror plane symmetry is a way to generate 
a strong topological insulator. 

FIRST PRINCIPLE CALCULATIONS 

We performed density functional theory (DFT) |[33l |34l 
calculations including Hubbard U and SOC using the 
full-potential linearised augmented-plane-wave (FP-LAPW) 
method as implemented in the elk code 11351 . The local density 
approximation (LDA) for the exchange and correlation energy 
functional in the Ceperley-Alder |[36l form parametrized by 
Perdew and Zunger 071 was employed. We used the "around 
mean field" (AMF) scheme [38| for the double-counting- 
correction. That is, to correct our DFT+U calculations for 
the Coulomb repulsion already present in the DFT Hamilto- 
nian. We treated up to 3d of Sr, up to 5 s of Ir excluding 4f, 
and Is of O with the radial Dirac equation, while the scaler 
relativistic approximation is used to include the SOC for the 
higher states in the second variational step 1591 . 
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FIG. 2: Some representative band structure diagrams of orthorhom- 
bic perovskite oxides for (a) magnetic metal (MM) at U = 1.5 eV 
and a = 0, (b) semimetal (SM) at U = 1.5 eV and a = 1.5, and (c) 
magnetic insulator at U = 2.5 eV and a = 1.5. The bands near the 
Fermi energy are denoted by red color, and remains knotted near U 
in SM phase. 



To confirm that our main results are robust to the choice of 
double counting correction, we have also computed the band 
structures using the "fully localized limit" (FLL) correction 
B0H421 near the phase boundary. We found that at large SOC, 
U c is essentially the same. However, for small SOC, U c is 



shifted towards a lower value than that found with AMF cor- 
rection in such a way that our main conclusion (the larger the 
SOC, the larger the U c ) does not alter. The phase boundaries 
obtained by these two different corrections are denoted by dif- 
ferent colours in the phase diagram shown in Fig. [3] and will 
be discussed below. 

To obtain the phase diagram of SrlrOs in the i/-SOC phase 
space, we tune the SOC term for the 5d orbitals of Ir atoms. 
Since the strength of the SOC increases sharply with the 
atomic number Z ( as Z 4 ) , it is much stronger for Ir (Z = 77) 
as compared to Sr (Z = 38) or O (Z = 8). This means that Ir 
contributes almost exclusively to the SOC energy in Srlr03. 
This allows us to safely tune the SOC for all valence states, 
because its effect on Sr and O atoms does not count much. A 
scaling factor a in the SOC term of the Hamiltonian is intro- 
duced in the second variational step ||39l . This way, we can 
enhance the effect of SOC by taking a > 1 or reduce it by 
taking a < 1 . For instance, a = would mean no SOC at all, 
while a = 1 is the atomic SOC in Ir atoms. A small magnetic 
field is used to set the quantisation direction for the angular 
momentum. This field reduces exponentially to zero during 
the self consistency iterations so it has no other effects. 

In the FP-LAPW method, the real space is divided into 
spheres around the atoms (muffin-tins) and interstitials else- 
where. In the present calculations, the muffin-tin radii 1.86 
a.u. , 2.08 a.u. and 1.51 a.u. are used for Strontium (Sr), Irid- 
ium (Ir) and Oxygen (O), respectively. The basis set consists 
of APW functions with angular momentum / up to 8 and plane 
waves with cut-off energy equal to 23 1 .3 eV. The number of 
empty states in the basis set in the second variational step was 
10. The Brillouin Zone integrations were performed using a 
3x3x3 grid, which is equivalent to using 10 points in the 
irreducible part of the Brillouin Zone. This works well, given 
that the primitive unit cell of orthorhombic perovskite Srlr03 
is almost four times bigger than that of the cubic structure with 
only one formula unit. We checked the k-grid convergence in 
the metallic phase using 8x8x8 grid. We only used U for 5d 
orbitals of Iridium. 



BAND STRUCTURES AND PHASE DIAGRAM 

The octahedral crystal field splits the bands derived from 
the of-orbitals of transition metal atoms into high energy e g 
and low energy t%g groups. Due to the distortion of octahedra, 
there are twelve tjg and eight e g bands (each band is doubly 
degenerate due to time reversal symmetry). Fig. refrep-bands 
shows band structures for various values of U and SOC de- 
noted in the inset. The crystal field gap between e g and tig is 
evident for all cases, and only bottom of e g bands are shown 
in the plots. 

When a — that corresponds to the absence of SOC, a 
ferromagnetic order is present, and t2 g bands are all mixed 
as shown in panel (a). In contrast, when a - 1.5, panel (b) 
and (c), the ti g bands form two groups, the higher four half 
filled bands originate from the J e ff = 1/2 denoted by red 
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colour ( the lower two bands near F point are mainly J e ff = 
3/2 though), and the lower eight completely filled bands from 
the J e ff = 3/2. Increasing the SOC increases the splitting 
between the J e ff =1/2 and J e ff = 3/2 bands. 




FIG. 3: The phase diagram of orthorhombic perovskite oxides in 
the (/-SOC plane. Three phases for U up to 4 eV and a = - 2 are 
Magnetic Metal (MM), non-magnetic Metal or Semimetal(M/SM), 
and Magnetic Insulator (MI). The colour circles show the points for 
which calculations have been performed and magenta, green and red 
denote MM, M/SM, and MI, respectively. Small magenta circles 
are FLL results showing MM phase where AFM gives M/SM phase. 
The solid line separates two phases connected via a first order phase 
transition (where blue line is obtained by AMF while magenta line is 
by FLL), whereas the dotted line is the phase boundary for a second 
order phase transition. 



For smaller U of panel (b), the phase is nonmagnetic 
semimetal, where four J e ff = 1/2 bands are near the Fermi 
level forming small pockets of Fermi surface. While in non- 
magnetic semimetal (SM) and magnetic metal (MM) phases, 
panel (a) and (b), there is a finite density of states at Fermi 
energy, the band topologies are very different in these two 
phases. In the non-magnetic semimetallic (SM) phase, the 
bands at the fermi energy cross near U point resulting in a line 
node, and the magnetisation is zero everywhere in this phase. 
In the MM phase, as well as in the magnetic insulator (MI) 
phase shown in panel (c), there is no such band crossing. In 
both these phases, Ir atoms have finite magnetic moments with 
a long range order. Increasing U, keeping the same strength 
of SOC, leads to a metal-insulator transition at a critical U c , 
where the insulating state such as panel (c) has an interest- 
ing magnetic ordering pattern. Since the time reversal sym- 
metry is broken due to the magnetic ordering, there are eight 
J e ff =1/2 bands in this phase as displayed in panel (c). A 
further discussion about the magnetic ordering pattern will be 
presented below. 

These three electronic phases shown in Fig. [2] are found in 
the £/-SOC phase diagram; (?) M/SM, (it) MM, and (Hi) MI. 
The overall phase diagram in U vs. SOC is presented in Fig. [3] 
where M/SM, MM, and MI are shown in green, magenta, and 



red. M/SM is connected to MM and MI via a first order phase 
transition whereby the magnetisation jumps from zero to a fi- 
nite value along with a sudden change in the band structure 
topology. On the other hand, MM and MI transform into one 
another continuously with opening or closing up of a band 
gap- 
Let us discuss the phase diagram by checking along dif- 
ferent cuts. First vertical cuts, i.e., changing U for a given 
a. When a = 0, the system remains a pure ferromagnetic 
metal at all U. This results from a large density of states at the 
Fermi level leading to a Stoner ferromagnet. Tuning SOC to 
finite but still small values (for a < 0.3), U interaction does 
not make any difference, and system stays in the magnetic 
metal phase even for very high values of U (for U close to 
5eV, it becomes ferromagnetic insulator, which is not shown 
here). However, as the SOC does not favour a pure ferro- 
magnetic ordering, it turns the magnetic ordering pattern to a 
slightly non-coplanar order with a large ferromagnetic com- 
ponent. In contrast, for a > 0.3, increasing U induces a first 
order phase transition from non-magnetic semimetal to mag- 
netic phases. Whether the magnetic phase is metal or insu- 
lator depends on the strengths of both U and a. The phase 
boundary separating the non-magnetic semimetal phase from 
the two other phases, MM and MI, is shown as a solid line in 
Fig. [3] For a > 1.1, increasing U transforms M/SM directly 
to MI, while for 0.3 < a < 1.1, increasing U changes the 
phase from non-magnetic metal, to magnetic metal followed 
by magnetic insulator. 

Let us explore the phase diagram using horizontal cuts - 
changing a for a fixed U. For small U, increasing a leads 
to a first order phase transition from magnetic metal to non- 
magnetic metal/semimetal phase. The critical value of a, a c , 
at which this transition takes place stays between 0.2 - 0.4 for 
< U < 1 .5 eV. This is rather expected, as SOC disfavours 
spin density wave ordering within a weak coupling theory. 
Fig. [3] shows results for a up to 2. As can be seen, there is 
no further phase transitions by increasing a. We checked this 
for a up to 5. 

For U > 1.5 eV, a c increases sharply with U with an in- 
creasing separation between the bands at the Fermi level in 
magnetic metal phase. It is also interesting to note that for 
2 < U < 2.35 eV, the system undergoes a change in phase by 
increasing a from magnetic metal to magnetic insulator, and 
then into non-magnetic semimetal, i.e, a re-entrance of metal- 
licity (metal-insulator-metal by changing SOC for a given U). 
For U > 2.35 eV, increasing a transforms magnetic metal 
smoothly to magnetic insulating phase with opening up of a 
band gap. The higher the value of U, the lower is the value of 
a for this transition. The phase boundary between these two 
phases is shown as a dotted line in Fig. [3] 

As can be seen in Fig. [3] the dotted line have the expected 
phase boundary curvature in the U-SOC plane. That is, the 
critical value of U for MM to MI transition decreases with in- 
creasing a since the phase space of magnetic insulator should 
be wider as U increases. However, the most unexpected re- 
sult presented in Fig. [3] is the opposite phase boundary cur- 
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FIG. 4: The J eff = 1/2 bands of Srlr0 3 close to U in the XURS- 
plane computed with LDA (without +U) at a = 1 (atomic SOC). 
These bands form two pairs of touching cones shown in yellow and 
brown. The two pairs are interpenetrated into each other forming 
a circular nodal line at the fermi energy. When the time reversal 
symmetry is broken (as the magnetic ordering occurs), these band 
crossings disappear and a band gap form. 
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FIG. 5: The bandwidth W of top two bands at fermi level as a 
function of a for various values of U. W decreases with a and U in 
MM/MI phase (U > 2.5 eV) as expected but increases in SM phase 
(U< 1.5 eV). 



vature for transition to magnetic phases from non-magnetic 
metallic phase, the solid line. It shows that a stronger SOC 
requires a stronger electron-electron interaction to transform 
non-magnetic metal to magnetic insulator. The origin of this 
unexpected behaviour is likely to be the special band topology 
in the semimetal phase as described below. 

Fig.[4]shows the band structure near U in the XURS-plane. 
The four J e ff = 1/2 bands form two interpenetrated pairs of 
cones, each pair consisting of a lower and a higher band (a 
yellow and a brown). One pair touches below the fermi level 
while the other above it, forming two Dirac-like points and 
a circular line node in the XURS-plane at the fermi energy. 
Due to this node, there is an extremely small density of states 
near the Fermi level, which in turn requires a high Hubbard 
U to splits these cones resulting in a magnetic insulator. We 
propose that this is the main mechanism whereby Srlr03 is a 
semimetal with a small carrier density, different from its sister 
compounds Sr 2 Ir0 4 dESrrE) and Sr 3 Ir 2 7 l|46ti48l. 

We also checked the bandwidth of the upper two bands of 
J e ff = 1/2 at the Fermi level when SOC is large enough to 
separate J e tt = 1/2 top two bands from the rest (except at 
the nodal points). W is plotted in Fig. [5] against a for various 
values of U for a > 0.5, since a < 0.5, J e /f = 1/2 is not well 
defined. Contrary to the expectation ll20ll . the bandwidth in- 
creases with a and U in the non-magnetic semi-metallic phase 
most likely due to a steeper slope of Dirac node that confirms 
our conclusion above. Whereas in the magnetic phases (the 
plots for U = 2.5 - 5 eV), W decreases with a and U as 
expected. This makes us believe that the transition from the 
non-magnetic to magnetic phases in S1T1O3 is controlled by 
the electronic state of the non-magnetic semimetal where the 
bandwidth is not relevant. This is further supported by the 
fact that the semimetal phase has a special band topology as 
described below. 



J. M. Carter et al showed in Ref. 1 32 ] using a tight binding 
model that this line node is proteced by the lattice symmetry. 
In other words, any term that opens up a gap near U-point 
should break either time reversal, inversion, or Pbnm lattice 
symmetry. For example, it was shown that when the mirror 
symmetry between the two layers of Ir02 in Srlr03 is bro- 
ken with a staggered potential, this line node changes to a 
point node. When the strength of this staggered potential is 
increased beyond a critical value that takes the node to the 
R point, a change in the topology of the bands occurs owing 
to the inversion of the J e ff = 1/2 bands at R. The system is 
turned into a strong topological insulator when this happens. 
Further increase of this staggered potential leads to inversion 
of the bands at Z point, changing the band topology back to 
trivial, making the system a band insulator. 

MAGNETIC ORDERING PATTERNS 

As discussed above, when a = 0, there is a pure ferromag- 
netic (FM) order in the MM phase at all values of U. This 
happens down to U = 0, because the electron-electron repul- 
sion is not completely absent even at U = 0, partly due to its 
imperfect removal in the 5a! orbitals of Ir and partly due to 
the presence of many other occupied states in the system. The 
magnitude of the magnetic moment of Ir depends on U. It in- 
creases with U from 0.38yu B at U = to 0.95yu B at U = 4 eV, 
where it is almost saturated — a rather expected behaviour. 
A small contribution to the magnetisation also comes from O 
when a ~ 0. 

As we move away from a = 0, system develops a canted 
antiferromagnetic (CAFM) order. At smaller U, transition 
from FM to CAFM is more gradual leaving a net ferromag- 
netic component. This behaviour persists up to a ~ 1 . For 
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FIG. 6: The magnetic structure of orthorhombic perovskite oxides 
at U = 2 eV and a = 0.2 (a), U = 4 eV and a = 0.2 (b), U = 4 
eV and a = 1 (c). In (a), the system has a canted antiferromagnetic 
order with a large ferromagnetic component. In (b), the system has 
a canted antiferromagnetic order with negligible ferromagnetic com- 
ponent. The average magnetisation per unit cell is also zero in (c), 
but the moments are aligned in a very different way. 



higher a a very small ferromagnetic component develops in 
the magnetic insulator phase at higher U. This is expected 
since an antiferromagnetic order in the insulating phase low- 
ers the energy via virtual hopping of electrons to the nearest 
neighbour with the oppositely aligned spin. A small ferromag- 
netic component is then due to an effective Dzyalonshinsky- 
Moriya interaction as found in Sr 2 Ir0 4 HH [T7] @9]|. Fig.[6] 
shows the magnetic structure at a = 0.2 as U is changed from 
2 to 4 eV, and at U — 4 eV as a is changed from 0.2 to 1 . The 
quantisation axis is set along the ;t-axis. 

For any two nearest neighbour Ir atoms in different layers 
along c (i.e, those with yellow and green or red and blue ar- 
rows in Fig.|6]l the components of moments along the y and 
the z-axis are always cancelled out. The size and orienta- 
tion/direction of individual moments depend on the values of 
a and U as does their sum or the total moment per unit cell. 
As can be seen in Fig.[6ja), at U — 2 eV and a = 0.2, the mo- 
ments are almost co-planer with a large ferromagnetic com- 
ponent. Fig.[6jb) shows the magnetic order at U = 4 eV and 
a = 0.2. It is clear from this figure that raising the strength of 
Coulomb interaction at finite a suppresses the ferromagnetic 
component. The magnetic order at the same value of U (4 eV) 
and a higher a, a — 1, is shown in Fig.|6jc), where a stronger 
SOC has changed the orientations and reduced the sizes of the 
individual moments. 



DISCUSSION AND SUMMARY 

The SOC is an essential ingredient in numerous exciting 
phenomena including spintronics and topological insulators. 
However, in transition metal oxides with 3d orbitals such as 
high temperature cuprates, the SOC has been ignored, while 
the strong correlation represented by Hubbard interaction de- 
termines their physical properties. Very recently, iridates with 
5d-orbital has been a topic of much attractive research due 
to intriguing combined effects of the SOC and Hubbard in- 



teraction. It was found that the SOC in iridates is unusually 
strong, which differs from other 5d compounds such as Re- 
oxides [50] even though the atomic SOC should be similar for 
Ir and Re. 

A set of iridates is the perovskite iridates forming Ruddles- 
den Popper series from single layer Sr2lr04 to three dimen- 
sional structure SrlrOs. While both single layer and bilayer 
iridates exhibit a magnetic insulating behavior, SrlrO^ shows 
a metallic phase with a small number of charge carriers. Given 
that the SOC and Hubbard interaction are local, their strengths 
should be similar in this series, and thus it was suggested that 
the bandwidth should control metal-insulator transition as the 
number of layers changes in perovskite iridates. [20 1 Here we 
show that the metallicity is innate to the lattice structure of 
three dimensional orthorhombic perovskites in addition to a 
large SOC. Due to this combined effect, there are tiny hole 
and electron Fermi pockets with small density of states, which 
in turn makes Hubbard interaction less efficient in SrlrO^. 
Due to strong SOC, the magnetic field dependence of phys- 
ical properties would be interesting to study. 

We investigate an overall phase diagram of the orthorhom- 
bic perovskite structure (space group Pbnm) in U vs. SOC 
using density functional theory. The computation is based 
on Srlr03, where tuning U and SOC (by changing a) al- 
lows us to explore other possible phases nearby non-magnetic 
semimetal in isostructural systems. Three phases - non- 
magnetic metal/semimetal, magnetic metal, and magnetic in- 
sulator - were found by tuning U and SOC. At smaller a, a 
magnetic metal is always found, which is similar to SrRuO^. 
While Ru 4+ has 4 electrons at the outer shell and thus the 
chemical potential is different from SrRh03, the bands near 
the Fermi level are strongly mixed leading to a similar phe- 
nomena. Indeed, earlier electronic calculation on SrRu03 re- 
ported it a ferromagnetic metal. At a > 0.3 and U < 1.5eV, 
the system becomes non-magnetic metal which resembles the 
ground state of SrRh03. Indeed, our computations of the elec- 
tronic structure of SrRh03 shows that it is similar to the one 
found at a = 0.4 close to the instability towards magnetic 
metallic phase. While Rh and Ru are next to each other in 
the periodic table, our results imply that the SOC must have 
a stronger effect on SrRh03 than SrRu03, and agree with an 
earlier suggestion that SrRh03 is near a magnetic critical point 
|@H6]. Increasing a further, the bands near the Fermi level 
changes to semimetallic-like, and a stronger U c is required 
for a magnetic insulator. The shape of phase boundary be- 
tween the non-magnetic semimetal and the magnetic insulator 
is emerged from a line of Dirac node leading a small density 
of states near the Fermi level. A tight binding approach for a 
series of Sr„ + ilr„03„ + i has found the same conclusion that U c 
is larger for n — oo than n — 1 or n — 2. ABU 

In summary, we have studied the interplay between the 
SOC and Hubbard interaction in orthorhombic perovskite ox- 
ide with the point group symmetry of Pbnm. Three different 
phases were identified. A magnetic metal with a finite fer- 
romagnetic component found in smaller SOC at all values of 
U investigated in this study. Increasing the SOC leads to a 
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phase transition to a non-magnetic metal for small U, and to 
a magnetic insulator for large U. The detailed band structures 
near the Fermi level in these phases strongly depend on the 
strength of the SOC rather than U, unless the interaction U 
leads to another magnetic phase. Our study may be useful in 
understanding different ground states found among isostruc- 
tural perovskites including SrRuC>3, SrRh03, and Srlr03. It 
also provides a microscopic mechanism for semimetallic be- 
haviour in Srlr03 distinct from its sister compounds, Sr2lr04 
and Sr 3 Ir 2 07. 
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